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We determine the nonperturbative anisotropic parameter of the gauge action in the quenched approximation 
with less than 1% accuracy using the Sommer scale measured by the Liischer-Weisz algorithm or smearing tech- 
nique. We also study the nonperturbative 0(a)-improvement of the quark action. The bare quark anisotropy 
is determined using the masses from the temporal and spatial directions. For the determination of the 0{a) 
improvement coefficients, we apply the Schrodinger functional method. 



1. Introduction 

Anisotropic lattices whose temporal lattice 
spacing Or is finer than the spatial one have 
become a powerful tool in various subjects of 
lattice QCD simulations. Among these appli- 
cations, computations of heavy-light matrix ele- 
ments |1I2I3I4| require the most accurate param- 
eter tuning, which should be performed nonper- 
turbatively. In this paper, we report the status 
of our project to develop the anisotropic lattice 
framework for such precision computations with 
accuracy of a few percent level [S] . 

Here we briefly summarize our strategy. For 
precise computations of heavy-light matrix ele- 
ments, we need a framework of the heavy quark 
in which one should be able to (i) take the con- 
tinuum limit, (ii) compute the parameters in the 
action and the operators nonperturbatively, (iii) 
and compute the matrix elements with a mod- 
est computational cost. The anisotropic lattice 
is a candidate of such framework, if a method 
which fulfills the above condition (ii) is provided. 
Our expectation is that on anisotropic lattices 
the mass dependence of the parameters becomes 
so mild that one can adopt coefficients deter- 
mined nonperturbatively at massless limit. We 
also need to control all the systematic errors in 
the continuum extrapolations. Feasibility studies 
performed so far for the level of 0(10%) compu- 
tations are encouraging for further development 
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|3I4| . We therefore investigate calibration proce- 
dures at the accuracy less than one percent, both 
for the gauge and quark actions in the quenched 
approximation at ^ = 4 [S]. 

2. Calibration of gauge field 

To achieve a few percent accuracy in the fi- 
nal results, the tuning of parameters must be 
performed at much less than this accuracy. As 
the goal of present work, we intend to determine 
the anisotropy parameters at 0(0.2%) level. The 
elaborated work by Klassen the 0(1%) level 
calibration for the Wilson action, is therefore no 
longer meets the present condition. For more 
precise calibration of the gauge field, we need 
to measure the static quark potential very accu- 
rately. For this purpose, we adopt the Liischer- 
Weisz noise reduction technique ;F as well as the 
standard smearing technique while applied in the 
anisotropic plane. The former method can dras- 
tically reduce the statistical errors while requires 
larger memory resources than the latter. 

We define the renormalized anisotropy 
through the hadronic radii tq measured in the 
coarse and fine directions. Since we carry out 
the continuum extrapolation in terms of the lat- 
tice scale set by rp , the renormalized anisotropy is 
kept fixed during the extrapolation. This avoids 
the systematic uncertainties due to the anisotropy 
which may remain in the continuum limit. 

Figure ^ shows a result of calibration at /3 = 
5.75. The top panel shows the result for the 
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Figure 1. Determination of the gauge field 
anisotropy at (3 = 5.75. 

static potential determined with the Liischer- 
Weisz technique at 7g = 3.072. The renormalized 
anisotropy is determined with 0.2% accuracy. 
A linear fit of the results at several values of 7g 
determines 7^ for which = 4 holds with 0.2% 
accuracy as displayed in the bottom panel. 

In Figure 12 the results at several values of /? 
are collected. Although the precisions of the re- 
sult with potential determined with the standard 
smearing technique are still not enough, sufficient 
precision is achieved at /3 < 6 where we used 
the L-W method. Improvement of calculation at 
/? > 6 and global fit analysis are in progress. 

3. Calibration of quark field 

Our heavy quark formulation basically follows 
the Fermilab approach but is formulated on 
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Figure 2. Result of calibration of gauge field. 

the anisotropic lattices |1I8| . The quark action is 
represented as 

Sf = J2^ix)Kix,yMy), (1) 

x,y 

K{x,y) = Sa:,y - Kr (1 - j4:)U4:{x)S ^^g ^ 

+ (l + 74)C/i(a;-4)<5,_4^^' 
+ ir + %)uHx-i)6^^-^y 

-KaCE^i(yAiFAi{x)5x,y 

-rKcrCBj2i>j'^ijPiji^)^x,y, (2) 

where Kq- and Kt are the spatial and temporal 
hopping parameters, r the spatial Wilson pa- 
rameter and ce and cb the clover coefficients. 
For a given k^, in principle, the four parameters 
7i? = Kt/k^, r, Ce and cb should be tuned. We 
can set r = 1 without loss of generality |1I7| . 

We must calibrate 7^, ce, and cb to the level 
which enables computations of matrix elements 
within a few percent accuracy. We also need to 
perform the nonperturbative renormalization of 
the operators such as the heavy-light axial cur- 
rent. The nonperturbative renormalization tech- 
nique jlO| is one of the most powerful methods to 
perform such a program. Following our strategy, 
this technique can also be applied with a little 
modification for the anisotropic lattice. 

We perform the calibration of 7i?, c^, cb, and 
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Figure 3. Tuning of at /3 = 6.1. 
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Figure 4. Determination of jp a.t P = 5.75. 

the renormalization coefficients of the axial cur- 
rent along the following steps. (1) Tuning of c^; 
by Schrodinger functional method. (2) Calibra- 
tion of 7i? (and if possible) by requiring the 
physical isotropy conditions for mps and my in 
the coarse and fine directions on lattices with T,L 
~ 2 fm. (3) Determination of cs, Hc and the 
renormalization coefficients of the axial current 
by Schrodinger functional method. We also need 
to verify that the systematic errors are under con- 
trol by calculating the hadron spectra and the 
dispersion relations and by taking the continuum 
limit. It is also necessary to verify that the tuned 
parameters in the massless limit is also available 
in the heavy quark mass region. 

To verify the feasibility of the step (1), we de- 



termine ce with fixed values of 7f and cb- At 
P = 9.5, the method is successfully applicable, 
and the result for cb is close to the 1-loop mean- 
field value. Figure O shows the result at /3 = 6.1 
(flo- ~ 2 GeV). The tuned value of ce is obtained 
as the value at which AM, the difference of the 
quark mass defined through the axial Ward iden- 
tity under different kinematical conditions, van- 
ishes up to 0{a^) effects. The result of c^; is larger 
than the tadpole improved tree level value. It is 
also found that AM is not sensitive to the change 
of Cb- 

Figure 01 shows the result for the step (2) at 
/3 = 5.75 on a 12^ X 24 X 96 lattice. The val- 
ues of is determined from the meson masses in 
the fine and coarse directions. We note that the 
result for is consistent with that from the dis- 
persion relation At this stage, the precision of 
jp is still not sufficient, while several techniques 
to reduce the statistical noise are yet to be tested. 
For the determination of c^, the ratio of the hy- 
perfine splittings in the fine and coarse directions 
is not feasible because of large statistical noise. 
Other procedures, such as Schrodinger functional 
method with boundaries in the coarse direction, 
are under investigation. 
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